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We develop a theory to model the degenerate two-photon state generated by the 50/50 Sagnac-loop source. We 

...... 

start with an interaction Hamiltonian that is capable of describing the interaction among the four optical fields 

°: 

' (non-degenerate pump and degenerate signal/idler), which reads: 

*° : 

t£ : Hj{t) = <*J dViE^E^E^E^+H.c) (1) 

t> : 

where a is a material constant that is characteristic of the optical fiber being used. The non-degenerate pump field 
^ is taken to be two synchronous pulses (denoted by subscripts pi and p2), copolarized and co-propagating down the 

\o '. 

fiber axis (denoted as z direction here), with central frequencies Cl p i and Cl p2 and envelope shapes E v \ and E p2 . 



Mathematically, they are written as below 



E p V = /rfw^i^Oe^ 1 ^-™" 1 ^-^ 4 , (2) 
E$ = I ' dw p2 E p2 {uo p2 )e lk ^)^^ e -^ > ( 3 ) 
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where lo p i = Cl p x + v p and io p2 = Cl p2 + v' p are the frequency arguments for the two pump fields. v p (y' ) denotes the 

CT 1 - 

frequency component within pi's (p2's) spectrum that deviates from its central frequency fi p i (Cl p2 ) by that amount. 

•rH . 

The phase tags e~ %lPlZ and e ~ l i p 2z are induced by pi's and p2's self-phase modulation (SPM) respectively, and are 
included in a straightforward manner. Note that P\ and P 2 denote the peak powers of pi and p2, respectively. 
The degenerate signal/idler field, with a center frequency at fl, is quantized according to [l[: 



= / dLJ a A(w a )ale- i[kM *-^ t] , (4) 
E { r } = I duiA^ale-^^-^, (5) 



where aj (j — s,i) is the creation operator for the j mode with frequency uij and wave-vector magnitude k(uij) = 
u> s = Cl + v s and u>i = CI + Vi represent the frequency of signal and idler photon, respectively, where v s 



c 



and Vi are the deviations for each photon's frequency from their central frequency fl. A(loa) = —it 57 — r is a 

slowly varying function of frequency and may be taken outside the integral. Now the interaction Hamiltonian may be 
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FIG. 1: Schematic of the Hong-Ou-Mandel experiment with the 50/50 Sagnac-loop identical photon source. FPC, fibre 
polarization controller; BS, beam splitter; OBPF, optical bandpass filter; Dl, D2, photon-counting detectors. Inset shows the 
spectral diagram of the non-degenerate pump and degenerate signal/idler fields. Q p i, pump-I (pi) central frequency; Q P 2, 
pump-2 (p2) central frequency; Q, signal/idler central frequency; a p , pump bandwidth; crn, OBPF bandwidth; A/2, central 
frequency difference between P2 and signal/idler (or signal/idler and PI); E s , Ei, Ei and E2, electrical fields before and after 
the BS, see text for details. 



expressed as 



Hi(t) = A J dzjdujsj dutalal J duj pl J dw p2 e~ i ^ +p ^ z E pl {u pl ) E p2 (u p2 ) 
exp{i[k(uj p i) + k(uo p2 ) - k{u s ) - k{uoi)]z - i{uj pl + uj p2 -oj s — u>i)t} , 



(6) 



where L is the length of the fiber, and A is an overall constant consisting of a, A(u>j) and the effective cross-section 
of the fiber A e g. 

The state vector at the output of the fiber can be calculated by means of first-order perturbation theory, namely, 



|*)out = |0) + 4 / H z (t)dt\0), 



(7) 



where the second term is the two-photon state that we seek, which we denote simply by |\&). After taking into account 
the fact that 



(8) 



we arrive at the following expression of the two-photon state: 



I*) = ^ / dz I <U. I dcj t atal\0) I <hj pl e~^ +p ^- 



ih 

where we have utilized the 5-function to simplify the integral over uj p2 , which also reinforces the energy conservation 

among the four interacting optical fields. 

To further simplify our analysis, we refer to the inset "Spectral diagram" shown in Fig. [1] which clearly illustrates 

the various system parameters by their corresponding mathematical symbols. These parameters correspond to the 

experimental settings schematically depicted in the main part of Fig. [TJ The two pump pulses are both assumed to 

be Gaussian-shaped with equal amplitude and equal bandwidth, i.e., 

("-"pi) 2 
E pl {ui) = E p e "i , 

(»-",2)' 

E p2 (uj) = E p e ^ , (10) 

where P\ = P 2 = P p cx E^ a 2 are the peak powers of the two pump pulses, and a p denotes their common optical 
bandwidth. By using Taylor expansion at the frequency 17 for the various fc's, we obtain 



Afc = k{u> p i) + k(uj s + uji — Upi) — k(uj s ) — k{u)i) 



k"(Q) 



A \ 2 /A 



(11) 



where we keep the expansion series to second-order dispersion only, which proves to be sufficient in most cases. 
The Wpi-integral in Eq. ^ can be rewritten as 

$(v s ,i>i,z) = J duj p iE p i{uj p i) E p2 (uj s + u>i - u; p i)e lAkz 

= El J d Vp e~ ^r'' + ^ K W + (»-*)(*^H» «] t (12) 

where in the last step we have used Eqs. (fT0|) and (fTTj) . and f3 2 is a shorthand for fc"(0). We are then left with the 
length integral 

,o 

Q{v„Vi) = J dz<P(v s ,v l ,z)e~ 2l ~< p r z , (13) 
and the two-photon state in Eq. (J9j) is reorganized into 



I*) = / / <<,',(){,',. 1>, + i>,)\ll + i',) . (14) 



where \ui) is a one-photon Fock state populated with a single photon of frequency to. We remark that the function 
Q(v s , Vi) is completely analogous to a(uj +uj e ) $>(lu , ui e ) in Eq. (9) in Ref. Similarly, \Q(v s , Vi)\ 2 can be interpreted 
as the probability distribution of the two-photon state 1\. However, the apparent symmetry of $(^ s , ^i, z), and thus 
Q(v s , Vi), with respect to its two frequency arguments results in qualitatively different behavior for the two-photon 
- & o m «„ t in Ref. Q, M > — c b * _ _, 

Further evaluation of $(2/ s , fj, 2) is made possible by using the integral formula from Ref. |2j, which deals with 
Gaussian integrals with complex arguments, resulting in 



("s+»j) A exp 
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(15) 



where A = f2 P 2 — Qpi is the central frequency difference between the two pump fields. We have thus formally obtained 
the expression for the two-photon state, which is given by Eq. (fT4"| or its following alternative version: 

2ttA 
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duj s I <koiQ{u) a ,Wi)\w s )\ui) , 



(16) 



where Q(oj s ,uii) is equivalent to Q(ui s — — ft) given by Eq. (fT5| . 

After obtaining the two-photon state, we are now ready to analyze the experiment shown schematically in Fig.[T] A 
variable delay St is inserted in one photon's path, before the two identical photons are recombined at the 50/50 beam- 
splitter (BS). As shown in Fig. [TJ if we denote the electric-field operators before the BS as Eg (t) and E\ + \t + Sr), 
then the field operators after the BS are given by 
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Ei + \t)+tEl + \t + Sr) 
iEW(t)+El +) (t + Sr) 



(17) 



where the vector nature of the field operators are ignored since they all share the same polarization, and 



E { +\t) oc / dwa, < a. 1< (w a , i )e- i "'-« t e 



(18) 



are the electric-field operators before the BS that include the shape of the OBPF, which is assumed to be Gaussian 



here. The coincidence-count rate registered by detectors Dl and D2 is given by 

p oo />oo 

Rc{8t) = dh dt 2 Pi2(h,t 2 ,ST), 
Jq Jo 
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(19) 



where 



Pl 2 (h,t 2 ,6r) = ^\E[-\t 1 )E!f'(t 2 )E^>(t 2 )E\ + '(t 1 )\^ 
= |<0|^ +) (i 2 )^ +) (*i)|*>| 2 
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(20) 



is the probability per pulse for coincidence detection between the two detectors. 

Eqs. HH), (fTT)) . (I18p . and (|20p . when plugged into Eq. fH)l) . after a simple but lengthy calculation, yield 



R c {5t) oc / du> s / du}i\F(u a ,Wi)\ 



l-e 



F(u! s ,uji) = Q(uJs,oJi)exp 



2o 2 



(21) 
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where we have used the fact that F(u> s ,u)i) = F{uji,uj s ) in obtaining the above results[7|. Its alternative version, 
written in terms of the difference- frequency v s and Vi, reads 



R c (5t) oc J dv s J dvi\F(v s ,Vi)\ 2 
F{v s ,Vi) = Q(!/ s ,i/j)exp 



2o 2 



(22) 



It turns out that further simplification is possible for the above Gaussian-filter case (by using again the formula 
from Ref. [2|), which we explicitly write out as the following: 



G{z) 



r 




/ dzi / 


J-L 


J-L 




r el* 2 


exp 


4(l+/3 



(23) 



exp 



exp 



- arctan(/3 2 zCTp) 



exp 



4(l + /3 2 2 z2 (T 4 ) 
V2n 2 P 2 a 2 cx p{i arctan 



i— 2i"/P p z 



§2 (zi-z 2 ) Op 

2 



(24) 
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(25) 



4 + p 2 ( Zl -z 2 ) 2 al 4 + p 2 ( Zl -z 2 ) 2 al 
One can, for instance, generalize the above result to investigate the effect that the filter has on the Hong-Ou-Mandel 
(HOM) dip. Due to its experimental relevance, we shall write out explicitly the formula that describes the case when 
the previously Gaussian-shaped OBPFs are replaced with two identical super-Gaussian filters, which reads 



R c (St) cx / dzi I dz 2 I dv s I dv l G{z 1 )G*{z 2 )e 2 "'i e 



CT p J-L 

exp 
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(26) 
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FIG. 2: Theoretical predictions vs. experimental results. Pink filled circles, experimental data; red solid curve, least-square 
fitting for the data; purple dotted curve, theory fitting from Eq. (|23[) for Gaussian OBPFs; blue dot-dashed curve, theory fitting 
from Eq. (|26[) for super-Gaussian OBPFs. Realistic values for the experimental parameters used in generating these curves are: 
L = 300 m, p 2 = -0.116 ps 2 /km, 7 = 1.8 x 10 -3 W _1 m~~\ P p = 0.36 W, c = 3 x 10 s m/s, A p i = 1555.92 nm, X p2 = 1545.95 nm, 
pump FWHM = 0.8 nm, signal/idler FWHM = 0.8 nm. 



where G{z) is given by Eq. 

We then generate two curves, from both Eq. (|23|) and Eq. (|26j) . to fit the experimental data obtained in the main 
text. The fitting results are shown in Fig. [21 As can be seen from the figure, both curves (blue dot-dashed and 
purple dotted) fit the experimental data remarkably well. The super-Gaussian fit appears to give a wider dip width 
(~ 8.0 ps FWHM) as compared to the Gaussian fit's result (FWHM dip width ~ 6.4 ps), which is commensurate with 
the fact that a super-Gaussian filter is spectrally narrower (and thus temporally wider) than its Gaussian counterpart 
with the same FWHM. A least-square Gaussian fit to the data (the red solid curve in Fig. (2), generated by the 
data-processing program, suggests that the HOM-dip visibility is ~ 94.3%. It has a FWHM dip width of about 7.2 ps, 
which is right in between the previous two fitting values. This may be explained by the fact that the real OBPF 
employed in the experiment is constructed by a cascade of a Gaussian filter and a super-Gaussian one, making its 
transmission spectrum somewhere in between. In contrast, the two theoretical fits agree on the ideally attainable dip 

l| that, as long as the 

two-photon probability distribution function[8j] is symmetric with respect to its two frequency arguments, the HOM 
dip can achieve a maximum visibility of 100%. 

There are many reasons that could explain the missing 5.7% visibility. For example, one might wonder whether 



visibility of 100%. This result also coincides with the theoretical understanding from Ref. 

.a 
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FIG. 3: Schematic drawing for investigating one of the possible scenarios for the less-than-unity HOM-dip visibility: spatial 
mode mismatch. See text for details. 

higher-order dispersion plays a role. But a straightforward Taylor expansion of the various fc's at the central frequency 

f2 [similar to Eq. to higher-order terms dismisses this hypothesis. In fact, the Ak quantity is always symmetric 

with respect to v s and Vi in our two-photon state due to the isotropic nature of the fiber. However, a little mismatch 

between the two OBPFs' spectrum will result in asymmetry of two arguments in F(v s , i/j), and will certainly cause a 

degradation of the HOM-dip visibility. Other candidates include: (i) The real-life BS's performance deviates from an 

2 R T 

ideal 50/50 BS, i.e., R + T = 1 and R^T . This gives rises to a corrective factor of ^ 2 ^ 2 *° ^ ne dip visibility 4]. 
When put in the measured values (R = 0.474, T = 0.526), it gives a 99.4% corrective coefficient, (ii) There might be 
some remaining ^2002 component due to the non-ideal alignment of the 50/50 Sagnac loop, which leads to degradation 
of the dip visibility 6]. (iii) The existence of other unsuppressed noise photons, such as Raman photons and single- 
pump FWM photons, could also degrade the attainable dip visibility, (iv) The spatial modes of the two photons are 
not exactly matched at the BS. A simple calculation, as we will carry out explicitly below, shows that a small angular 
mismatch between the two photons' paths distinguishes between the coincidence-generating amplitudes (TT and RR). 
As a result, they are not completely cancelled after the BS, which correspond to the remaining coincidences at the 
HOM dip. An angular mismatch as small as 3 x 10 _5 rad is required to bring the dip visibility down to 94.3%. The 
details of the calculation go as follows. Suppose at the BS the two photons intersect at a small angle 8, as depicted 
in Fig. [3j The two amplitudes can each be written in the Fourier-optic language as: 

4 

4 



V'tt = --e 12 ^ 2 



V'RR = ^e^Pi. (27) 

7T0T 
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Here d is the diameter of the lenses used to couple light into fibre, pi.2 represent the two-dimensional coordinates at 
the lens planes (perpendicular to the paper), and |/o| = |/o| = are the projected wave- vector magnitudes in the 
lens planes for the two off-axis waves. The overlap between the two amplitudes, which is proportional to the HOM-dip 
visibility [a], is given by 

where O(cc) (it obtains the value 1 when \x\ < 1 and everywhere else zero) represents the effective areas of the 
lenses, and Ji(x) is the first-order Bessel function. When we put in realistic values for d — 5mm and A = 1.55 /im, 
and demand that Eq. (128]) yields 0.943, we obtain a numerical solution for 9 ~ sin 9 ~ 30 /irad. From the above 
calculation, we can thus see that spatial mode mismatching has the highest likelihood of contributing to the missing 
5.7% HOM-dip visibility. 



Ji(7rrf|sin6>|/A) 
7rd|sin6»|/(2A) 



(28) 
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